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Abstract: Recently several works have appeared in the literature that addressed the problem of 
Freeze Out in energetic heavy ion reaction and aimed for a description based on the Boltzmann 
Transport Equation (BTE). In this paper we develop a dynamical Freeze-Out description, starting 
from the BTE, pointing out the basic limitations of the BTE approach, and the points where the 
BTE approach should be modified. 



PACS numbers: 25.75.Ld, 25.75.-q, 24.10.Nz 

The Freeze Out (FO) is an important phase of dynam- 
ical reactions. The connection of the kinetic description 
of this process and the Boltzmann Transport liquation 
(BTE) raised considerable attention recently 0, Q 

The phase-space (PS) distribution of post FO particles 
can be obtained from kinetic FO calculations. These were 
performed earlier in one-dimensional models |^ Ij, M; |^ , 
where the dynamics was governed by two constants: re- 
thermalization and FO characteristic scales. 

The FO is a kinetic process and one would think it 
can be handled perfectly by using the BTE, which may 
describe both equilibrium and non-equilibrium processes 
in a 4-dimensional space-time (ST) volume element like 
FO layer. This work and [l| follows this approach. Our 
aim is to analyse the situation, discuss the applicability 
of BTE, and point out the physical causes, which limit 
the applicability of the BTE for describing FO. 

FO is usually assumed to happen on sharp 3- 
dimensional ST hypersurfaces. However, the FO-fronts 
or FO-layers are not necessarily narrow, but they have 
a characteristic direction (or normal, d?a^), and it is 
more realistic to assume a continuous, 4-volume FO in 
a layer (or domain) of the ST. At the inside boundary of 
this layer there are only interacting particles, while at the 
outside boundary hypersurface all particles are frozen out 
and no interacting particles remain. Then the ST volume 
element, d'^x, in the layer of interest can be converted into 
dl^x — *■ ds d^a^, where ds is the length element in the 
direction of the 4- vector d^a^ (= da^ for short), which 
can be time-like or space-like. Let us assume that the 
boundaries of this layer are approximately parallel, and 
thus thickness of the layer does not vary much. 

The emission or FO probability may depend on phys- 
ical processes, cross-sections, transition rates, and the 
actual PS distributions, f{x,p). Furthermore, / must 
be determined self- consistently during the FO, i.e. all 
conservation laws should be satisfied across this tran- 
sition layer and overall entropy should not decrease! 



iSifff,?,"^. 

a. Modified Boltzmann Transport Equation: We can 
derive the BTE from the conservation of charges in a 
ST domain, A^x, assuming the standard conditions [lOJ: 
(i) only binary collisions are considered, (ii) we assume 
"molecular chaos", i.e. that the number of binary colli- 
sions at position x is proportional to f{x,pi) x f{x,p2), 
and that (iii) f{x,p) is a smoothly varying function com- 
pared to the mean free path (m.f.p.). We have to take 
into account that particles can scatter into the PS volume 
element around p, or can scatter out from this volume el- 
ement, described by Gain- and Loss- collision terms in 
the BTE. In these terms we consider elementary binary 
collisions where in the initial state two particles collide 
with momenta pi and p2 into a final state of two particles 
with momenta p^ and p^. In case of the Gain term the 
particle described by the BTE, with momentum p (with- 
out an index), is one from the two final state particles, 
while in case of the loss term this particle is one of the 
initial state particles. This is indicated by the indexes 
of the invariant transition rate ^^. We integrate over 
the momenta of the other three particles participating 
in this binary collision. We use the notation for the PS 



integrals: u'Ds 



d-'pi d-'p2 d-'ps 
Pi P2 Ps. 



Shortening the notation 



further, by dropping the momentum arguments of the PS 
distributions and keeping the indices only, yields: 



P'^dj - i 1 ,2VJlf^Wf^ -\j 2V,J^fiW^ 
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p2 



(1) 



Now, aiming for the description of FO let us split up the 
distribution function, / = p + f^ , where f^ is the PS 
distribution function of the " free" or frozen out particles, 
while /* is the interacting component. 1, 3] As the parti- 
cles belonging to the free component may not collide any 
more, they do not appear in the initial state components 
of collision integrals! 

The gain term, fl /g Wfj , populates both the in- 
teracting, /\ and free, /■* , components. So, we intro- 



duce a FO probability, which 'feeds' the free component, 
Vf = V^'^{x,p), and the rest, (l—Vf), feeds the interact- 
ing one. Now, we can separate the two components into 
two equations. The sum of these two equations returns 
the complete BTE above: 






\ j 12V JifiVfWl^ , (2) 

-\j 12V JlflVfWl^ (3) 

\ J 12V jinwi^ ^\j 2v,,rf^w^^ . 



The free component does not have a loss term, because 
particles in the free component cannot collide, and so, the 
free component cannot loose particles due to collisions. 
The first term of the second equation is a drain term, de- 
scribing the escape or FO of particles from the interact- 
ing component. It is the inverse of the gain term for the 
free component. The last two terms are influencing the 
interacting component, and do not include the FO prob- 
ability factors! Thus, these two terms drive the interact- 
ing component towards re-thermalization. For these the 
relaxation time approximation can be used [J, |5|, Ifil • 

P''d,f - -i 1 12V JlfiVfWl^ + P° ^^^^ . (4) 

^ J Trel 

We see that the structure of the kinetic equations, used 
earlier |j,|j,IMI^, and the separation of the "escape" and 
"re-thermalization" terms come out in a simple straight- 
forward way from the BTE. 

However, the usual structure of the collision terms in 
the BTE is not adequate for describing rapid FO, in a 
layer which is comparable to the m.f.p. If we assume 
the existence of such a layer this immediately contra- 
dicts assumption (iii): the change is not negligible in the 
direction of da'^ . The assumption of "molecular chaos" 
is also violated in a FO process because number of col- 
lisions is not proportional with f{x,pi) x f{x,p2), but 
it is delocalized in the normal direction with f{xi,pi) x 
f{x2,P2)- (The fact that the FO is a delocalized kinetic 
process, was already used in ref. jlj when integrals along 
the path of propagating particles were introduced, but 
the consequences regarding the details of the collision 
terms and the validity of the molecular chaos assump- 
tion were not discussed.) 

Based on the above considerations, one might con- 
clude that the changes of the distribution function are 
mediated by the transfer of particles, and consequently 
only slowly propagating changes are possible. If the FO 
layer propagates slowly, thus, its normal, da^, must al- 
ways be space-like. This was a common misconception, 
where all " superluminous" shock, detonation, deflagra- 
tion fronts or discontinuities were considered unphysical 
based on early studies [ll|. However, it was shown re- 
cently, that discontinuous changes may happen simulta- 
neously in spatially neighbouring points, i.e. the normal 
of the discontinuity- hypersurface can be time- like |l2l |. 



This applies to the FO process also. Thus, the direction 
of characteristic or dominant change, da'^, may be both 
space-like and time-like in the FO process. 

From the all processes mentioned above (i.e. shocks, 
detonations, deflagrations etc.) the FO is the most spe- 
cial one, because the number of interacting particles is 
constantly decreasing as the FO proceeds, correspond- 
ingly the m.f.p. is increasing. In fact, it reaches infin- 
ity when the complete FO is finished. This means that 
we can not make the FO in finite layer of any thickness 
smooth enough to be modeled with the BTE. It is also 
obvious that if FO has some characteristic length scale, 
it is not proportional with the m.f.p. only, because the 
m.f.p. increases as the density of interacting component 
becomes smaller, while the FO becomes faster in this 
limit, so its characteristic scale should decrease. 

To describe that the PS distributions change rapidly 
along the FO direction, wc can introduce the Modified 
Boltzmann Transport Equation (MBTE): 



p^d^ff{x,p) 
P^d^fix.p) 



\ j 12ViV fWl^ n{Xi,pi)n{x2,P2) , 

-\ j 12V^V fWl^ f{{xuPl)r2{.X2,P2) 
^0 -f Gq J 



'^rel 



(5) 



where Xk is the origin of colliding particles, i.e., the ST 
point where the colliding particles were colliding last, 
Xk ^ X - UkTk, Tk is the collision time, u^ = {"fk,JkVk), 
7 = l/Vl — v^ and Vk = Pk/Pk- This is an essential mod- 
ification if the PS distribution has a large gradient in the 
ST. This gradient defines a ST 4- vector, characterizing 
the direction of the process du'^ . In ref. |3] the direction 
dcr'^ is also introduced, however, it is not discussed why 
and it is not connected to the delocalization of the BTE. 

The symmetries and the assumption of local molecular 
chaos leads to the consequence that local conservation 
laws can be derived from the original BTE, i.e. d^/T^^'^ ~ 
and dfj^N'^'- = were T and N are given as momentum- 
integrals over the single particle PS distribution. After 
we delocalize the equations, this leads to new type of 
transport terms. The FO process results in additional 
particle and energy currents in the FO direction. 

In the simplest case of infinitely narrow FO layer (FO 
hypersurface) between in and out domains the distribu- 
tion function is given by / = finQin + fout^out and com- 
mon BTE can be used everywhere except for the FO 
hypersurface. This case was actually studied without in- 
troducing MBTE in ref. U and the additional currents 
discussed above led to the possibility of the new type of 
Phase Transition - the so called three flux discontinuity, 
see 2j for more details. 

The molecular chaos is also essential in the Boltzmann 
H-theorem. If we relax this condition, the solution of 
MBTE does not converge to a stationary thermal distri- 
bution, which enables us a simple thermo- and hydrody- 
namical treatment of the physical problem. 



b. Approximate Kinetic Freeze Out Models: Let us 
present the schematic derivation of the kinetic FO model 
used by some of the authors ear her. This represents only 
one particular possibility and the general MBTE equa- 
tion can be solved or approximated in other ways also. 

If the ST distribution is non-uniform and the direction 
of steepest gradient can be clearly identified, one may re- 
place one (or more) of the integrals over (Ppi (or (fip2) by 
ST integrals over the origins of the incoming particle (s), 
d'^x', requiring that the particle reaches the ST point, x, 
when needed. This requirement determines p'^ for a given 
x''^. It is reasonable to assume that after converting some 
of the integrals to ST integrals and performing them, we 
get an effective FO term reflecting the properties of local 
PS distribution, transition rate, and ST configuration. 




FIG. 1: The plot of one of the last collisions at x, shown 
in the spatial cross section of the FO layer. Particles arrive 
from positions xi and X2 to point x with momenta pi and 
P2- Within the FO layer of thickness, L, the density of in- 
teracting particles gradually decreases (indicated by shading) 
and disappears at the outside boundary, Si (thin line) of the 
layer. R.h.s. from this boundary there are no interacting 
particles. Particles can reach x from a region closer than the 
mean free path (m.f.p. indicated by the dashed line), but 
only from places where the interacting particle density is still 
not zero, i.e., mostly from the left. The inside boundary of 
the FO layer, S2 (thick line) indicates the points where the 
FO starts. Left of this line there is only interacting matter 
and the FO probability is assumed to be zero for collisions 
happening in the interacting region. 

Let us study the first equation of |SJl, without per- 
forming the integrals in a formal way, rather illustrat- 
ing the procedure giving a better insight to the prob- 
lem. When we are in the FO layer close to the bound- 
ary of complete FO, we have to calculate here the 
collision rate . According to the MBTE this depends on 
the PS distribution of the incoming particles at their ori- 
gins, fl{xi,pi) X fl{x2TP2)- Assume that the FO points 
in the direction x (to the right), then on the r.h.s. of the 
collision point density of particles is low or zero, while 
on the l.h.s. it is larger, close to the pre FO value (see 



Figure^. Most probably both incoming particles ar- 
rive from the left, where still there are many interacting 
particles and very few or none from the right. Conse- 
quently most outgoing particles leave to the right. Thus, 
the collision rate at x depends on the conditions around 
X — A, i.e. deeper inside the interactive matter. So, the 
collision rate is higher than the conditions at x could se- 
cure! It does not go to zero even if we are at the outside 
boundary of the FO layer. Si, because particles still can 
arrive from the left. 

Let us execute two of the PS integrals for one incoming 
and one outgoing particle, / — p- — S^ 



Pl 



o^djf 



1 



iiVJUlVfWi 



pi 



Q2Vi fd'^p 



!^fl{xi,pi)VfWll , 



p\ 



(6) 



where Q2 = J — r 12(^2, P2) and V4 are invariant scalars. 
Eq. © resembles cq. (3.27) in ref. [13, but one of the 
incoming particle distributions, /2(3^2,P2) is integrated 
out, and leads to an integral quantity, Q2- This can 
be approximated by the invariant scalar density at X2, 
i.e., Q2 ~ ri2{x2) ■ Here V4 is not given explicitly, but 
can in principle be calculated based on the distributions, 
fi{xi,pi) and f2{x2,P2), and the transition rate, W. The 
resulting transition rate will then be averaged over par- 
ticles 2 and 4: Wf;^ . So, finally we obtain: 



p^d^ff{x,p) 



Q2Vi fd^pi 



Pi 



fl{xi,pi)VfW^^ 



(7) 



where xi = {ti,xi) = (i - 71^™/;., x - viJiTcou.) = 
X — ui TcoiL, (with the notations Vi — Pi/pi and 
""1 = (71 J 71^1)) indicates that the particle arrives at x 
starting earlier from another point xi. As the first rough 
approximation we can take that TcooI. is a mean collision 
time. The more general way to take into account particles 
coming from the different distances with corresponding 
propability factors is presented in |l5|. In a rapid dy- 
namical process the distributions at x and xi are not the 
same, as it was discussed above. Now, eq. {7} can be 
integrated either in the xi -space, or the pi -space. 

In addition, the FO probability, Vf, may include inte- 
grated information about the FO process, e.g. the prob- 
ability not-to-coUide with anything on the way out. This 
should depend on the integral number of interacting par- 
ticles on the way out. 

For simplicity let us assume small angle scatterings, 
and the propagation of a single particle W^^ « w^ 5{p — 



Pi) 



then 



p^djfix,p) 



Q2V4 



r{xi,p) Vf wt 



(8) 



where P{xi,p) can be expressed in terms of f^(x,p) and 
the 4- velocity, u'', of the particle with momentum p. 



Let us now consider the FO situation, where we have 
a directed process in a layer. The dominant change hap- 
pens in the direction of the normal of the FO hyper- 
surface, da^^ {da^ da^ = ±1) and it is negligible along 
the hypersurface of the front. This actually means that 
p^dufix.p) « {p^dcTf^) da^diy f{x,p). Inserting the above 
equation into ||Sll yields to a kinetic equation describing 
the directional derivative of the distribution function in 
the direction of the dominant change, da'^ , as 



da^d,ff{x,p)^nx^,p)p:,,, 



(9) 



where f''{xi,p) reflects the system properties one m. f. p. 
earlier than the collision point and the escape probability 
Pesc = P*sc{^: Pi Pi d.'^i ^1 ^/)- The derivation above did 
neglect several details and features, however, reflects the 
basic structure of ad hoc kinetic FO models H, Hi IM Hi - 

c. Escape Probability: Here we just present briefly a 
direct estimate for the escape probability [Ijl, detailed 
introduction and analysis of which will be given in [Ij] . 

The escape probability includes the collision rate and 
FO probability. The probability not-to-coUide with any- 
thing on the way out, reasonably should depend on the 
number particles, which are in the way of a particle mov- 
ing outwards in the direction p/p, across an FO layer of 
estimated thickness L (representing the fact that we have 
finite number of particles on the way out to collide with 
PI). If we are in this FO layer in a position x'^, there 
is still 7T^ distance ahead of us, where O is an an- 

cos ' 

gle between the normal vector and p/p. We assume then 
that the FO probability is inversely proportional to some 
power of this quantity 13]. Thus 



esc \ / - \ 

A(a;i) 



L 



L — xf^dan 



(cose)"e(p''da^), (10) 



where the power a is influencing the FO proflle across the 
front, and the cut factor is eliminating negative contribu- 
tions to FO. In papers [1,|^|MIS authors used a — 1, and 
modeled FO in an infinite layer {L — > oo), and a constant 
characteristic length A were used instead of A(a;i): 



P*^^ = e(p^dcr^)cose/A. 



(11) 



Comparing eqs. l(Tn|l and ((TTll one can see that we 
replaced the constant characteristic length A, which 
was clearly oversimplifying the situation, with A'(x) — 
A(xi)-^^ , (for daf^ — (0,1,0,0)) which contains two 



factors. The first is the collision rate, which is propor- 
tional with jyt^ ~ {n{xi)a) , and this does not tend 
to zero even if we reach the outside boundary of the FO 
layer, as this parameter is characteristic to the interior 
region at xi. The other is the generalized FO probability, 
which depends on the direction of the outgoing particle 
and on the number of interacting particles left in the way 
to collide with, i.e. oc -p^- Now the new characteristic 
length A' (x) gradually decreases as FO proceeds and the 
number of interacting particles decreases, and goes to 0, 
when the FO is finished, as it was discussed above. 

The simple angular factor, cos 8, maximizes the FO 
probability of those particles, which propagate in the 
direction closest to the normal of the layer, da^. The 
quantities, cosO = p^/|p| for FO in x-direction and 
cos = 1 for FO in t-direction ;6.] , are not Lorentz in- 
variant. Therefore to make our description completely 



invariant we shall generalize it to 






cosO 9]. 



So, we write the invariant escape probability, within 
the FO layer covering both the timelikc and space-like 
parts of the layer [13 , as 



P* 



L 



X{xi) \L — x>^d<7fj 






Qip^'da, 



(12) 

In refs. 0, p| the post FO distribution was evaluated 
for space-like gradual FO in a kinetic model with the es- 
cape probability (|ll|l . The post FO distribution depends 
on the details of the escape probability and on the level 
of re-equilibration of the interacting component. Bugaev 
assumed earlier |7| that the post FO distribution is a 
(sharply) " Cut- Juttner" distribution, but the above men- 
tioned model shows that this can only be obtained if re- 
equilibration is not taking place. The kinetic modelpro- 
vided an asymmetric but smooth PS distribution |j, |5[, 
while the escape probability (|12|l yields a somewhat dif- 
ferent, but also smooth PS distribution [13^ . 

d. Conclusions: We have shown that the basic as- 
sumptions of BTE are not satisfied in the FO process, 
and so, the equation should be modified if used for the 
description of FO. We have also shown that earlier ad 
hoc kinetic models of the FO process can be obtained 
from our approach in a fully covariant way, and FO in 
space-like and time-like directions can be handled on the 
same covariant footing. 
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